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Due to the Heisenberg uncertainty principle, various classical systems differing only on the scale 
smaller than Planck's cell correspond to the same quantum system. This fact is used to find a 
unique semiclassical representation without the Van Vleck determinant, applicable to a large class of 
correlation functions expressible as quantum fidelity. As in the Feynman path integral formulation of 
quantum mechanics, all contributing trajectories have the same amplitude: that is why it is denoted 
the "dephasing representation." By relating the present approach to the problem of existence of 
true trajectories near numerically-computed chaotic trajectories, the approximation is made rigorous 
for any system in which the shadowing theorem holds. Numerical implementation only requires 
computing actions along the unperturbed trajectories and not finding the shadowing trajectories. 
While semiclassical linear-response theory was used before in quasi-integrable and chaotic systems, 
here its validity is justified in the most generic, mixed systems. Dephasing representation appears 
to be a rare practical method to calculate quantum correlation functions in nonuniversal regimes in 
many-dimensional systems where exact quantum calculations are impossible. 
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The method that is described in this paper is based on 
two observations: the first is the fact that the relation- 
ship between classical and quantum dynamics is many to 
one; the second is the idea of shadowing of a perturbed 
trajectory by a nearby unperturbed trajectory. We start 
by explaining these two ingredients in more detail. 

Classical vs quantum dynamics. In the semiclassical 
(SC) approximation, quantum wave function is associ- 
ated with a classical surface (Lagrangian manifold) in 
phase space. SC evolution of the wave function is per- 
formed by classically evolving this surface and computing 
actions along the trajectories of the points of the surface. 
At the end, the surface is projected onto an appropri- 
ate coordinate plane. If we slightly distort the initial 
surface, individual trajectories will change exponentially 
fast. However, if the distortion is small enough, the orig- 
inal and distorted initial surfaces semiclassically corre- 
spond to the same quantum wave functions (their over- 
lap is ~ 1). Due to the unitarity of quantum evolution, 
the overlap of the two wave functions associated with the 
two evolved surfaces will remain ss 1 for all times. 

Shadowing. Because of the exponential sensitivity to 
initial conditions and because of the finite precision of a 
computer, computer generated trajectories in chaotic sys- 
tems are accurate only for a logarithmically short time. 
As a result, it was not clear whether it makes sense at all 
to do computer simulations for longer times and whether, 
e.g., the fractal patterns seen in these simulations are 
real. The solution was offered by Hammel, Yorke, and 
Grebogi Q with the concept of shadowing which was 
later, in various settings, promoted to a theorem [HQ 



Their finding is that while a computed trajectory diverges 
exponentially from the true trajectory with the same ini- 
tial conditions, there exists an errorless trajectory with a 
slightly different initial condition that stays near ( "shad- 
ows") the computed one. 

While shadowing is extremely useful if the perturba- 
tion is random, which is the case for the roundoff com- 
puter errors, it works also if the perturbation is determin- 
istic, e.g., given by a precise change of the Hamiltonian. 
This is a completely different physical problem: Even 
if a computer could evolve trajectories exactly without 
any roundoff errors, we might be interested in compar- 
ing the trajectories with the same initial conditions, but 
of two slightly different Hamiltonians. It turns out that 
in chaotic systems these trajectories will again exponen- 
tially diverge. Below we rigorously show that the latter 
problem can be transformed to the former. 

Comparison of the dynamics of two slightly different 
Hamiltonians is the subject of a very rich and recently 
much studied subject of the sensitivity of classical and 
quantum motion to perturbations |4J, yj la, \A 1<J m ■ This 
sensitivity is best described by so-called fidelity M (t) : 
Classical (quantum) fidelity is the classical (quantum) 
overlap at time t of two initially identical phase-space 
distributions (quantum states \ip}), that were evolved by 
slightly different Hamiltonians H° and H e = H° + eV 
where e controls the strength of the perturbation. In 
Dirac notation. 



M(t) = \0(t)\' 



ifj e 



,+i(H°+eV)t/h -iH°t/h 
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Studying fidelity is useful in itself but also because 
many other, seemingly unrelated quantum correlation 
functions take the same mathematical form. For exam- 
ple, M(t) is equivalent to the Loschmidt echo — the sur- 
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vival probability of a state evolved first by H° for time 
t, then by —H e for time t. If we were to evaluate this 
overlap directly semiclassically 0, IE j we would have 
to find exact trajectories of H° and of H e , corresponding 
phases and prefactors, and then calculate the overlap. 
This would be extremely difficult in chaotic systems be- 
cause of the exponential growth of errors, exponentially 
proliferating contributing trajectories, and exponentially 
growing number of singularities. In fact, full SC treat- 
ment, i.e., evaluating the overlap itself by the stationary 
phase approximation, would yield anything but the right 
answer since we would be adding up a huge number of 
infinite singularities. Compared to the direct approach, 
the dephasing representation (DR) is simple, as can be 
seen from the following alternative expressions: 

DR {t) = {2irhy d J d d p' exp [iAS (r', p', t) /h] , (1) 
DR (t) = J dVexp [iAS(r',p',t)/H}\i>(r')\ 2 , 



DR (t) = J d d r' J dVexp 



-AS(r',p',t) 



Pw (r',p') , 



where the first expression is for position eigenstates |r'), 
the second for wavepackets with mean momentum p', 
and the third (and most general) for Wigner distribu- 
tions pw ■ Since Wigner distribution can represent mixed 
states, the third expression can describe purity fidelity 
and decoherence. Analogs of the first two expressions 
exist in momentum and other representations, and all 
can be derived from the Wigner form. In all expressions, 
the action difference AS is just the integral of the per- 
turbation along the unperturbed trajectory with initial 
conditions r', p', 

AS (r', p', t) = S e - 5° = -e f drV [r (t) , r] . 
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Roughly speaking, the motivation for DR is the following: 
instead of using the trajectory of H e with the same initial 
condition, we use the shadowing trajectory: a trajectory 
of H e with a slightly different initial condition, which 
remains close to the trajectory of H° up to time t. [That 
this is possible follows from the shadowing theorem (ST) 
as we show below.] All we have to compute is the phase 
difference along the trajectory. We avoid the singularities 
since the prefactors precisely cancel and we also need 
to evaluate a much smaller number of trajectories than 
in the standard SC treatment. We explain these claims 
below when we derive the approximation. 

First let us show that the problem of sensitivity of dy- 
namics to the change of Hamiltonian is equivalent to the 
problem of shadowing of numerically noisy trajectories. 
For simplicity we show this for a general two-dimensional 
symplectic map, which is nothing else but a discretiza- 
tion of continuous-time dynamics of a Hamiltonian sys- 
tem. We need three standard definitions that hold for 
any (not-necessarily symplectic) map / : x„ i— > x n+ i. 



Definition. A true trajectory {x„} n=a of a map / 
satisfies x Jl+ i = /(x n ) for a < n < b. 

Definition. {x ra }^ =Q is an e-pseudotrajectory of a map 
/ if |x n+ i - /(x„)| < e for a < n < b. 

Definition. The true trajectory {x„}^ =a S-shadows 
{x. n } n=a on a < n < b if |x„, — x n | < S for a < n < b. 

Anosov ^(j an d Bowen pj| showed that for uniformly 
hyperbolic maps, noisy trajectories can be shadowed for 
an infinitely long time by true trajectories. Most physical 
systems, however, are not uniformly hyperbolic, which 
led to a series of investigations for more general systems. 
By now, there exist theorems for both continuous differ- 
ential equations and discrete maps, and they require the 
satisfaction of various assumptions 0, 0, 0- The gen- 
eral conclusion is that shadowing in nonuniformly hyper- 
bolic systems only works for a finite time ts, which de- 
creases with increasing perturbation. Below we relate the 
problem of Hamiltonian perturbations to the problem of 
the shadowing of noisy trajectories for symplectic maps, 
without considering detailed assumptions of a particular 
ST. 

Let q n and p n be the position and momentum at time 
n and let x„ = (q n ,p n ). Let the continuous system be de- 
scribed by a Hamiltonian H i = H Q +eV = p 2 /2+W+eV, 
where W is the unperturbed potential and eV is the per- 
turbation. Further assume that V and W are differen- 
tiable and |VV(<7, f)| < 1 for all q. The corresponding 
symplectic map f : x„ h x n+ i is given by 

Pn+i = Pn ~ VW(fo,n) - eVV(q n ,n) , 

Qn+l = Qn + Pn+1- 

Lemma 1. The true trajectory {x„}^ =Q of f e is an 
e-pseudotrajectory of /°. 



/°(x„) 



/ e (x w )-/°(x n ) 



Proof |x n+ i 
W (q n )\<e. 

Lemma 2. Let the noise amplitude be S, 
the numerically computed trajectory {x„}^ =a of /' 



i.e. 



let 
be 



Then {x n } n=a is a (S + e) 



a 5-pseudotrajectory of f e 
pseudotrajectory of f°. 

Proof |x n+ i - /°(x n )| 

|x„+i - r(x„) + r (x„) - /°(x„)| < | Xn+1 - r( x „)i + 

|/ £ (x„)-/°(x„)| < S + |/ e (x„)-/°(x n )| < S + e, 
where we used Lemma 1 in the last inequality. 

Lemma 1 shows that Hamiltonian perturbations are 
equivalent to random noise as long as further assumptions 
of the particular ST are satisfied. Lemma 2 shows that 
we can combine Hamiltonian perturbations with random 
noise before applying the ST. 

The author believes that the DR expressions Q are 
the fundamental representations of fidelity because they 
provide (to his knowledge) not only the only accurate but 
simply the only way to calculate fidelity semiclassically 
for longer than the logarithmic time. DR can be thought 
of as the Van Vleck propagator in the mixed position- 
momentum representation, so that the Van Vleck deter- 
minant equals det((9pfi na i/<9p') = 1 and the Maslov in- 
dices cancel. However, since most authors are familiar 
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only with the Van Vleck propagator in the position rep- 
resentation, we show explicitly how DR can be derived 
by a careful treatment of this propagator. 

The Van Vleck propagator in position representation, 



K SC {v",v'-t) = ]T (2^)- d/2 cj /2 exp ( l -S 3 



is the standard SC approximation of the quantum prop- 
agator. The sum is over all classical trajectories connect- 
ing point r' at time t = to point r" at time t = t, Sj 
is the classical action, 

S j (r" ) r';t)= f drL [r (r) ,f (r) ,r] , 
Jo 

and Cj = | dct(d 2 Sj/dr"dr')\ is the Van Vleck determi- 
nant equal to the classical transition probability. There 
are three major problems with this approximation. First, 
the number of trajectories in a generic system grows ex- 
ponentially with time. Second, there is a growing num- 
ber of singularities in this expression at conjugate points. 
Third, direct evaluation of this expression would involve 
an expensive root search, i.e., finding for each trajectory j 
the initial momentum leading to the final position r". 
In improved SC methods such as the initial value repre- 
sentation, the root search and singularities are avoided, 
at a cost of replacing the sum over trajectories j by an 
integral over initial conditions 0, 0] . 

Forgetting for a moment all these problems (as have 
forgotten most authors discussing fidelity semiclassically 
0,H, with the exception of Ref. where fidelity is com- 
puted semiclassically for a logarithmically short time), let 
us use Kgc to express fidelity amplitude: 

0(t)« / drK* s * c (r,r';t)K sc (r,r';t) 



(The superscript tells us which Hamiltonian is used.) We 
could proceed in several ways: we could try evaluat- 
ing this integral by the stationary-phase method. That 
would give a singular result as explained in the text above 
Eq. (Q) because there can be many coalescing stationary 
phase points. (For disbelievers the author recommends 
trying it.) Or we could evaluate the integral numerically, 
which might smooth out some of the singularities in the 
two integrands by integration, but would be virtually im- 
possible in practice. What Jalabert and Pastawski p| 
suggested was using only the diagonal terms (j — j e ), 
justifying this by the claim that the off-diagonal terms 
cancel out in the average over initial states (or real- 
izations of the perturbation, since they calculated the 
average fidelity). There are several problems with this 
approach. Since H° and H e are different, after a long 
enough time it will not be possible to distinguish terms 
as diagonal and off-diagonal. Also it is not obvious that 
the off-diagonal terms should cancel since these are in 
majority. 



However, this separation is possible in DR. For the 
trajectories of -ff e , instead of using the precise initial con- 
ditions, we infinitesimally adjust the initial conditions so 
that the trajectories of H e shadow those of H°. While 
Kdr = ^sc remains unchanged in DR, 



K^ R {vy,t) = Y;(^h)- d /\c^ 



exp 



Note that because we are using the shadowing trajectory, 
the Maslov index and the Van Vleck determinant are the 
same as for K^ R , only the action is different (since this 
trajectory is a true trajectory of H e ). Now there is an 
exact one-to-one correspondence between the terms of 
the two propagators. Moreover, if we use the Loschmidt 
echo picture, only the diagonal terms form a continuous 
trajectory in phase space (momentum at time t is contin- 
uous), and so semiclassically, only these terms survive. It 
turns out, that in DR, neglecting the off-diagonal terms 
actually improves the approximation. This can be seen, 
e.g., for e = 0, where we get 0_D«(i) = 1 only if we ne- 
glect the off-diagonal terms (this was noted already by 
Jalabert and Pastawski i). As a result, we have 

DR (t) = J d d rJ2(2nhy d C°exp(iAS j /h). 

3 

Realizing that C° = | det(<9p'/<9r)|, we obtain DR Q for 
position eigenstates |r'}. 
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FIG. 1: Comparison of the exact quantum calculation (dots) 
and dephasing representation (dashed-dotted line) of fidelity 
M{i) for the standard map (n = 1000, 1000 classical trajec- 
tories used). Top: mixed phase space (k — 0.8, e = 5 x 10 -3 ). 
Bottom: chaotic phase space (k — 10, e = 2x 10~ 3 ). 

DR was first numerically tested in chaotic systems 
which are "closer" to being uniformly hyperbolic 0, 0] • 
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Motivated by the ST for nonhyperbolic systems here we 
also apply the method to mixed systems. Figure [T] shows 
fidelity decay of initial position states (q = 0.4) for the 
Chirikov standard map (W — kV = —kcosq, H = l/lim) 
in cases with mixed and chaotic dynamics. The success 
of DR is striking: using only 1000 trajectories in the 
chaotic case (bottom), the approximation is excellent at 
time t = 50 where there are ~ 10 20 SC contributions 
in the standard approach H, H, @ • DR works since for 
chaotic systems, only the statistics of actions matters, 
which can be reproduced by a smaller ensemble [|| Q- 
Counterintuitivcly, experts in the field of chaos might 
be more surprised by the results for mixed phase space 
(top). The reason is that in both integrable and chaotic 
systems there exist different simplifications. In mixed 
systems with both invariant tori and chaotic regions, nei- 
ther simplification works throughout phase space. DR 
still reproduces even the fine details of M(t) in Fig. ^ 

While there is a rigorous prescription how to verify the 
applicability of the ST, it is beyond the goal of this pa- 
per. It is done by numerically finding, for each trajectory, 
rigorous bounds on several dynamical quantities [3] ■ The 
goal of this paper was to reveal the surprising connection 
between classical shadowing and quantum mechanics and 
to present a generally applicable SC method. We there- 
fore omit here the rigorous check of the assumptions of 
the ST for the specific example in Fig. ^ A rough esti- 
mate ts ~ e -1 / 2 for the shadowing time from a conjecture 
in Refs. [3, HQ gives t s ~ 14 and 23 for the top and bot- 
tom parts of Fig. ^ respectively. This estimate should 
be used with caution since it does not even depend on k. 
Figure ^ suggests that most trajectories are shadowable 
well beyond this time, because otherwise it would be vir- 
tually impossible that a method based on the ST, using a 
large number of quite randomly interfering waves should 
mimic so accurately the exact quantum solution. 

A more thorough description of the universal numerical 
success of DR in the Gaussian, algebraic, Fermi-Golden- 
Rule, Lyapunov, and intermediate regimes, and for more 
general perturbations, will be presented elsewhere |14| . 



Furthermore, it can be shown that considering statistics 
of action differences in DR leads to a simple unified the- 
ory of the four regimes [hi, ITE) . It should be noted that 
before the work in Ref. |6j , due to the exponential prolif- 
eration of SC contributions, numerical evaluation of M(t) 
for longer times was usually only done using exact quan- 
tum propagation on a grid |8( ■ For the standard map that 
we used, the most successful SC calculations were done 
up to time t < 10 where there were fewer than 10 4 con- 
tributing classical trajectories In many-dimensional 
systems, DR seems to be the only method to evaluate 
fidelity because the complexity of exact quantum propa- 
gation grows exponentially with dimensionality. 

To conclude, DR is an efficient method that can de- 
scribe quantum-mechanical decays purely in terms of de- 
phasing. This strangely suggests that in some cases it 
is possible to describe the decay of classical overlaps by 
interference fl4T |. In particular, we have rigorously shown 
that DR is accurate for systems and times where ST 
applies. The many-to-one relationship between classi- 
cal and quantum dynamics on which DR is based has 
been exploited before 0, . While the replacement- 
manifold method [TtI ] is very simple and gives excellent 
results in a variety of problems, it requires finding the 
replacement manifolds. DR goes one step further: here 
the fact that a slightly distorted initial surface exists is 
sufficient. It is not necessary to find it explicitly. 

The author conjectures that DR should be applica- 
ble to other classes of correlation functions of the type 
{ip\e Ai e B \ip) where operator A differs only slightly from 
B, and A and B do not necessarily correspond to real- 
time quantum evolution. Possible applications may 
include temporal and thermal correlation functions in 
many-body systems in condensed matter and chemical 
physics, and in the general wave scattering in disordered 
media. 

The author would like to thank E.J. Heller and W.H. 
Miller for discussions and the Mathematical Sciences Re- 
search Institute for support. 



[1] S.M. Hammel, J. A. Yorke, and C. Grebogi, J. Complex. 
3, 136 (1987). 

[2] C. Grebogi et al, Phys. Rev. Lett. 65, 1527 (1990). 

[3] T. Sauer and J.A. Yorke, Nonlinearity 4, 961 (1991); S.- 

N. Chow and K.J. Palmer, J. Complex. 8, 398 (1992); 

B.A. Coomes, H. Kocak, and K.J. Palmer, J. Comput. 

Appl. Math. 52, 35 (1994); T. Sauer, C. Grebogi, and J. 

A. Yorke, Phys. Rev. Lett. 79, 59 (1997); H. Kantz et 

al., Phys. Rev. E 65, 026209 (2002); W. Hayes and K.R. 

Jackson, SIAM (Soc. Ind. Appl. Math.) J. Numer. Anal. 

41, 1948 (2003). 
[4] A. Peres, Phys. Rev. A 30, 1610 (1984). 
[5] R.A. Jalabert and H.M. Pastawski, Phys. Rev. Lett. 86, 

2490 (2001). 

[6] J. Vam'cek and E.J. Heller, Physical Review E, 68, 
056208 (2003). 



[7] J. Vamcek, Ph.D. thesis, Harvard University, 2003, 
http: //physics .harvard. edu/Thesespdf s/ vanicek.pdf 

[8] P. Jacquod, PG. Silvestrov, and C.W.J. Beenakker, 
Phys. Rev. E 64, 055203 (2001); 

P. Jacquod, I. Adagideli, and C.W.J. Beenakker, Euro- 
phys. Lett. 61, 729 (2003); 

F. M. Cucchietti, D.A.R. Dalvit, J. P. Paz, and W.H. 
Zurek, Phys. Rev. Lett. 91, 210403 (2003); 

G. Benenti, G. Casati, and G. Veble, Phys. Rev. E 67, 
055202(R) (2003); 

G. Benenti, G. Casati, and G. Veble, Phys. Rev. E 68, 
036212 (2003); 

B. Eckhardt, J. Phys. A 36, 371 (2003); 
T. Kottos and D. Cohen, Europhys. Lett. 61, 431 (2003); 
P.G Silvestrov, J. Tworzydlo, and C.W.J. Beenakker, 
Phys. Rev. E 67, 025204(R) (2003); 



5 



M. Hiller, T. Kottos, D. Cohen, and T. Geisel, Phys. Rev. 
Lett. 92, 010402 (2004); 

T. Gorin, T. Prosen, and T. H. Seligman, New J. Phys. 
6, 20 (2004); 

T. Prosen and M. Znidaric, New J. of Phys. 5, 10 9 (2003); 
T. Prosen and M. Znidaric, quant -ph/0401142 (2004); 
D. V. Bevilaqua and E. J. Heller, |nli~ CD/0409007 
(2004). 

[9] N.R. Cerruti and S. Tomsovic, Phys. Rev. Lett. 88, 
054103 (2002). 

[10] D.V. Anosov, Proc. Steklov Inst. Math. 90, 1 (1967). 
[11] R. Bowen, J. Diff. Eqns. 18, 333 (1975). 



[12] W.H. Miller, J. Chem. Phys., 53, 3578 (1970). 

[13] W.H. Miller, J. Phys. Chem. A 105, 2942 (2001). 

[14] J. Vamcek, invited presentation, March meeting of Amer- 
ican Physical Society, Montreal (2004). 

[15] J. Vamcek, e-print: |http : //arXiv\ org (2004). 

[16] L. Kaplan, New J. Phys. 4, 90.1 (2002). 

[17] J. Vamcek and E.J. Heller, Phys. Rev. E 64, 026215 
(2001); J. Vamcek and E.J. Heller, Phys. Rev. E 67, 
016211 (2003); J. Vamcek and D. Cohen, J. Phys. A 36, 
9591 (2003). 



